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Abstract 

Two-phase flow in pneumatic two-fluid fuel noz- 
zles was experimentally investigated to determine 
the effect of atomizing-gas density and gas mass- 
flux on liquid-jet breakup in sonic-velocity gas- 
flow, Oropsize data were obtained for the following 
atomi zing-gases: nitrogen, argon, carbon dioxide 

and helium. They were selected to cover a gas 
molecular-weight range of 4 to 44. Atomizing-gas 
mass-flux ranged from 6 to 50 g/cm^-sec and four 
differently sized two-fluid fuel nozzles were used 
having orifice diameters that varied from 0,32 to 
0.56 cm. 

The ratio of liquid-jet diameter to SMD, D 0 /D32. 
was correlated with aerodynamic and liquid-surface 
forces based on the product of the Weber and Reynolds 
number. We Re, and gas- to- 1 1 qu i d density ratio, 

Pg tp\. To correlate spray dropsize with breakup 
forces produced by using different atomi zing-gases, a 
new molecular-scale dimensionless group was derived, 

p^ V^/ giig , and used to obtain the following 
expression: 

V°32 = 4.33xl0- n (we Re Pg/p,) 0 ' 4 ^, V^/ giJg ) 0 * ?5 
2 2 3 

where We Re = p^D^/op^ , p 1 is liquid viscosity, 

a is surface tension, V c is acoustic gas velocity, 
V m is RMS velocity of gas molecules, g is acceler- 
ation of gas molecules due to gravity, and pg is 
gas viscosity. 

From this expression for liquid-jet breakup 
in two-fluid atomizers, it is evident that 

°32 ~ (pgVq) 1 ’ 33 . This expression agrees well with 

atomization theory for liquid-jet breakup in high- 
velocity gasflow, i.e., in the regime of aero- 
dynamic stripping. Also, it was found that at the 
same gas mass-flux, pgV c , helium was considerably 
more effective than nitrogen in producing small- 
droplet sprays with SMD's in the order of 5 pm. 

Nomenclatur e 

b dropsize parameter in Nukiyama-Tanasawa 

expression, cm 

c dropsize parameter in Rosin-Rammler expres- 
sion, cm 

D q characteristic drop diameter measured for 
entire spray, cm 

D) diameter of ith drop, cm 

0 V ,5 volume median drop diameter, cm 

D 31 volume-linear mean drop diameter, 


D 32 Sauter mean drop diameter, 

HinD?/]Ti n 0 i’ cm 

k correlation coefficient for Eq. (1) 

k' correlation coefficient for Eq. (3) 

N n exponent for Nukiyama-Tanasawa dropsize dis- 
tribution expression 

N r exponent for Rosin-Rammler dropsize distribu- 

tion expression 

n number of droplets 

R gas constant, 8.31x1 0 7 ergs 

Re Reynolds number, D 0 pgV c /p 

V fluid velocity, cm/sec 

v volume fraction of droplets with diameters 

less than or equal to x 

W weight flow of fluid, g/sec 

2 

We Weber number, D„p V to 

o r g c 

x droplet diameter in dropsize distribution 

expression, cm 

p absolute viscosity, g/cm sec 

p fluid density, g/cm3 

0 surface tension of liquid relative to air, 
dynes/cm 

Subscripts : 

c acoustic 

g gas 

1 liquid 

m molecular 

o orifice 

w water 

Introduction 

Liquid-fueled gas turbine and rocket combustor 
performance can be markedly Improved with two-fluid 
fuel nozzles capable of producing large numbers of 
small fuel droplets, i.e., fuel sprays that have a 
large liquid surface-area to volume ratio. This 
improvement in combustor performance Is due to the 
fact that large surface-area sprays yield high 
vaporization and burning rates. The very small 
droplets produced by this type of atomizer vaporize 
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fast enough so that it may also he used to form 
premi xed-prevaporized fuel-air mixtures suitable 
fo.r application in catalytic combustors. 

As reported in Ref. 1, an investigation was 
previously conducted using four differently sized 
two-fluid nozzles. It was found that reciprocal 

Sauter mean diameter, . could be correlated with 

the product of the Weber number, Reynolds number, 
and gas- to- 1 i qu i d density ratio raised to the 0.44 

power as follows: ~ (we Re p This 

-1 33 

relationship shows that ~ v c ’ where V c is 

acoustic gas velocity. The exponent agrees with 
that predicted by atomization theory for liquid-jet 
breakup in high velocity gasflow. 

Numerous investigators have studied the effect 
of fluid properties on liquid-jet atomization. Some 
of their results are given in Refs. 2 to 4. Of all 
the fluid properties that affect atomization, it was 
found in Ref. 1, that acoustic mass-flux, pgV c , had 
the greatest effect on the fineness of atomization 
obtained with two-fluid nozzles. However at pres- 
ent, the independent effect of gas density, p g , on 
spray dropsize has not been very well established in 
the literature as indicated in Refs. 5 and 6. There- 
fore, in the present study, the effect of gas density 
on spray dropsize was investigated. To accomplish 
this, acoustic gas-velocity was held constant while 
gas density was varied. In addition, a wide range 
of acoustic gas-velocities was investigated by using 
four different atomizing gases having molecular 
weights ranging from 4 to 44. The gases used were 
helium, nitrogen, argon and carbon dioxide. 

In the present study, two-phase flow in fuel 
nozzles was experimentally investigated to deter- 
mine the effect of acoustic gas-velocity and mass- 
flux on characteristic spray dropsizes produced by 
liquid-jet breakup in high velocity gasflow. Char- 
acteristic dropsizes were measured with a scattered- 
light scanner developed at the NASA Lewis Research 
Center. 7 The sprays were sampled at a distance of 
2.2 cm downstream of the fuel-nozzle orifice, to 
minimize the loss of small droplets due to vaporiza- 
tion and dispersion. The sprays were injected into 
a low velocity, 5 m/sec, stream to aid in transport- 
ing very small droplets through the laser beam. 
Small-droplet sprays with Sauter mean diameters 
varying from 3 to 30 pm were produced by the atomi- 
zation of small-diameter liquid jets in the regime 
of aerodynamic-stripping, i.e., high-velocity gas- 
flow. Liquid and gas pressures were varied from 
0.2 to 1 .0 MPa. 

Apparatus and Procedure 

A two-fluid atomizer was mounted horizontally 
in the open duct as shown in Fig. 1. Air supplied 
at ambient temperature, 293 K, passed through the 
24 cm inside diameter test section before being 
exhausted to the atmosphere . The test section was 
1 m in length and airflow rate in the test section 
was measured with a 5.08 cm diameter orifice. To 
aid in transporting small droplets through the 
laser beam, an airstream velocity of 5 m/sec was 
maintained in the test section. Figure 2 shows a 
schematic drawing of the spray and Instrumentation. 

Four pneumatic two-fluid atomizers were fabri- 
cated according to the diagram illustrated in Fig. 3 


with orifice diameters ranging from 0.32 to 0.56 cm. 
The atomizer was operated over pressure ranges of 
0.2 to 1.0 MPa for both water and the atomizing gas. 
Sprays were injected downstream into the airflow 
just upstream of the duct exit. The sprays were 
sampled at a distance of 2.2 cm downstream of the 
atomizer. 

Water at a temperature of 298 K, measured with 
an I.C. thermocouple, was axially injected into the 
airstream by gradually opening a control valve until 
the desired flow rate of 3.15 g/sec was obtained as 
indicated by a turbine flowmeter. The atomizing gas 
was then turned on and weight flowrate was measured 
with a 0.51 cm diameter sharp-edge orifice. After 
the air, atomizing gas and water flowrates were set, 
the Sauter mean drop diameter was measured with the 
scattered-light scanner to characterize each spray. 
Exponents for the Rosin-Rammler and Nukiyama-Tanasawa 
dropsize distribution expressions were also deter- 
mined using the scattered-light scanner. 

The optical system of the scattered-light scan- 
ner shown in Fig. 2 consists of a laser beam expander 
with a spatial filter, rotating scanning slit, and a 
detector. Using a 4.4 by 1.9 cm laser beam, the 
instrument measures scattered light as a function of 
scattering angle by repeatedly sweeping a variable- 
length slit in the focal plane of the collecting 
lens. The data obtained is scattered-light energy 
as a function of the scattering angle relative to 
the laser beam axis. This method of particle size 
measurement is similar to that given in Ref. 8. 

In Ref. 9, it was found that measurements of 
points on a plot of scattered-light energy normal- 
ized to the maximum energy and plotted against scat- 
tering angle can be used to determine the Sauter 
mean diameter, the volume median, the volume linear 
and the 75-percent volume diameters as illustrated 
in Fig. 4. The dispersion of the size distribution 
is also determined from the points on the plot shown 
in Fig. 4. Also, it should be noted that this method 
of determining characteristic drop diameters and dis- 
persion of dropsizes can be used independent of the 
particle size distribution function as discussed in 
Ref. 9. For a typical measurement, the scan is 
repeated 60 times per second to average out any 
temporal variations in the energy curve. 

Spray pattern effects were minimized by measur- 
ing characteristic drop diameters using the entire 
spray. Reproducibility using this technique was 
within ±5 percent. Calibration was accomplished 
with five sets of monosized polystyrene spheres hav- 
ing diameters of 8, 12, 25, 50, and 100 pm. Since 
the present sprays were sampled very close to the 
atomizer orifice, they contained a relatively high 
number-density of very small drops. As a result, 
the light-scattering measurements required correc- 
tion for multiple scattering as described in Ref. 10. 

Drop size measurements were also corrected to 
include Mie scattering theory when very small drop- 
let diameters, i.e., <10 pm, were measured. 

Experimental Results 

Four atomizing gases, i.e., nitrogen, argon, 
carbon dioxide and helium were used in two-fluid 
fuel nozzles to study the effect of atomizing-gas 
density on spray dropsize. To demonstrate this 
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effect, several dimensionless groups were Effect of Atomi zi nq-Gas Properties on Liquid-Jet 

investigated. Breakup 


Effect of Atomi zinq-Gas Mass-Flux on Spray Dropstze 

The reciprocal SMD, was plotted against 

atomizing-gas flowrate per unit area as shown in 
Eig. 5. From this plot, the following relationship 
is obtained for the four atomi zi ng-gases , helium, 
nitrogen, argon and carbon dioxide: 

D ;]- k ( y \>)' 33 » 

where k is a correlation coefficient that is a 
function of the gas used to atomize water jets. 
Equation (1) can be rewritten in terms of gas mass- 
flux as follows: 



Values of the correlation coefficient, k, for each 
atomizing gas are given in Table 1. 

The four atomizing qases show good agreement 
with atomization theory,'' which predicts recipro- 
cal SMD to be directly proportional to gas velocity 
raised to the 1.33 power. This exponent also agrees 
well with the results previously reported in Ref. 1 
and the value of 11.5 obtained for the coefficient 
k is in good agreement with that reported in Ref. 7 
for two-fluid atomizers using nitrogen gas. 


The primary objective of this investigation was 
to study liquid-jet breakup in two-fluid atomizers 
by using four different atomizing-gases and corre- 
late reciprocal Sauter mean diameter, with 

atomizing-gas properties. As shown in Fig. 6, the 
product of the Weber and Reynolds numbers and the 
density ratio do not give a single correlating 
expression. In order to accomplish the objective 
of this study, the normalized reciprocal Sauter 
mean diameter, D 0 /D 32 , produced by liquid-jet 
breakup with two-fluid atomizers is also assumed to 
be a function of the RMS velocity of the gas mole- 
cules, V m , the acceleration of gas molecules due to 
gravity, g, atomizing-gas viscosity, p g , and liquid 
density, pi . Using dimensional analysis, the fol- 
lowing expressions are obtained: 

°o /D 32 = f< VV»V9> (4> 

where = ^3RT g /M g . 12 By rewriting Eq. (4), 
there results: 

°o /D 32 * f <P,> a <V%> C V d (5> 

The preceding equation is then expressed in terms of 
the mass-length-time system (where T is time: M, 

mass; and L, length) to give: 


A comparison of k values clearly shows the 
marked Improvement that is obtained in fineness of 
atomization when a low-density gas such as helium is 
used as the atomizing gas for two-fluid fuel nozzles 
operating at acoustic gas velocities. The attain- 
ment of a large liquid-surface area per unit volume 
of liquid, i.e., large values of reciprocal SMD. 

D^. with the use of helium is attributed to the low 

molecular weight and therefore high acoustic-velocity 
of helium. On the other hand, the higher molecular 
weights and lower acoustic velocities of argon and 
carbon dioxide make them much less efficient as 
atomizing-gases for two-fluid fuel nozzles operating 
at acoustic gas velocities. 



so that for: 


M, 0 - a + d 
T, 0 > -b -2c -d 
L, 0 » -3a +b +c -d 
which may be written as: 

a - -d 


( 6 ) 


Correlation of Normalized Reciprocal SMD with Dimen- 
sionless Groups 

In the study of liquid-jet breakup in high- 
velocity nitrogen gasflow that is reported in 
Ref. 1, it was found that the reciprocal of the nor- 
malized SMD, (D 32 /D 0 ) -1 , could be directly corre- 
lated with the product of the Weber and Reynolds 
numbers, We Re, multiplied by the gas- to- 1 i qu i d den- 
sity ratio, pg/pi, as follows: 

D o /D 32 = k ' (We Re p g /p l )0 ’ 44 (3) 

2 2 3 

where We Re = D o p g V c /pjO and V c is the acoustic 

velocity of the gas. Atomization data obtained for 
the four atomizing gases are shown in Fig. 6. From 
the slope of the plots, a correlation coefficient 
k‘ was determined for each atomizing-gas and val- 
ues of k' are recorded in Table 1. 


b = -2c -d = -3d 
c - 3a -b +d 

Substitution of these values into Eq. (5) gives: 

D o /D 32 - f ( p l V m / %) _d < 7 

To obtain a general correlating expression for 
the four atomizing-gases, the correlation coeffi- 
cient k‘ is plotted against p t V 3 /gp g . From the 

slope of this plot shown in Fig. 7, the following 
relationship is obtained: 

k ' ~ ( p l V >g) < 8 



where the exponent -d, derived from dimensionless 
analysis is equal to 0.75. Since this exponent is 
fairly large, the variables p g , p i , V m , and g 
have a very important effect on the liquid-jet 
breakup process. Microgravity studies are needed 
in which g is varied to determine if the exponent 

0.75 is correct for g. 

To obtain a single correlating coefficient for 
the atomizing-gases , values of D 0 /D 32 are P'otted 
against the product of the dimensionless groups 
given in Eqs. (3) and (8), as shown in Fig. 8. Thus 
the following expression is obtained for the four 
atomizing gases: 


D 0 /D 32 - ^xlO-^He Re 

(9) 

Since We Re(pg/pj) = D o( p g V c) 3/p l p l CT ’ ' ' s ev 'dent 
that the reciprocal Sauter mean diameter is propor- 
tional to V^' 1 , p°' 31 , Pg 0 ’ 75 an d gas molecular 

weight raised to the -0.46 power. For the liquid 

-1 -0 44 -0 44 

properties, D,, is proportional to p, ' , o 

0 31 

and pl| . These exponents give fairly good agree- 
ment with those predicted by atomization theory, as 
shown in Table 2. Therefore Eq. (9), although it 
was derived for water-jet breakup, should give 
fairly good results in predicting fuel-jet breakup 
in Mach 1 gasflow with two-fluid nozzles. 

Correlation of Exponents for Dropsize Distribution 
Expressions 


The scattered-light scanner gave data for the 
exponent N r , which appears in the Rosl n-Rammler 
expression as follows: 33 


N x r N r 

dv _ r -(x/c) r (10) 

dx = cN f 

Oata were also obtained for the exponent N n , which 
appears in the Nukiyama-Tanasawa expression as 
follows : 3 3 


dv 

dx = T(6/N n > 


6/N N 

n 5 -bx n 
x e 


( 11 ) 


From a plot of the data obtained with the four 
atomizers, as shown in Fig. 9, the following rela- 
tion was determined: 

N = 2.8N 0-45 (12) 
r n 


which is the same as that derived in Ref. 1. Thus 
it was found that the relation between exponents N n 
and N r was not appreciably affected when the 
atomizing-gas molecular weight was varied from 4 
to 44. 


Summary of Results 

From the correlating expressions derived in 
this study, it is evident that: 


ORIGINAL F « -- 
OF POOR QoALJTY 


1. The effect of atomizing-gas density, p g , on 
spray dropsize is the same as that of acoustic 3 

velocity, V c , i.e., Dj 2 - ( Pg\/ C > 1 -33 - 

2. RMS gas molecular-velocity, V m , and acoustic 
mass-flux, pgV c , have the greatest effect of all 
gas- and liquid-phase properties on the process of 
liquid-jet breakup in high-velocity gasflow. 

3. The reciprocal SMD, produced with two- 

fluid atomizers was considerably greater for helium 
in comparison with nitrogen, at the same gas mass- 
flux. The superior performance of helium gas as 
compared with that obtained with nitrogen, is 
attributed to the relatively high acoustic and RMS 
molecular velocities of helium. 

Appendix - Characteristic Dropsize of Liquid Fuel- 
Sprays Simulated by Water-Jet Breakup 
in High-Velocity Gasflow 

Values of liquid properties given in Table 3 
were used in Eq. (9) to calculate reciprocal Sauter 
mean diameters for JP-5 and n-Octane fuel-sprays. 
These are plotted against gas mass-flux in Fig. 10. 
Two sets of plots are presented. The elevation of 
the upper set over the lower one shows the effect 
on spray surface-area that is obtained, according 
to Eq. (9), when helium gas is used in two-fluid 
atomizers instead of nitrogen gas. From these 
plots, values of k were calculated and are 
recorded in Table 3. 
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TABLE 1. - CORRELATION COEFFICIENTS FOR 
ATOMIZING GASES 


Correlation 

coefficients 

k, Eq. (1) 

k' . Eq. (3) 

86.5 

53 

11.5 

7.5 

6.3 

Tt.O 

6.3 

4.0 


LIQUID-JET PHYSICAL PROPERTIES AT 20 *C AND VALUES 
OF k FOR EQ. (1) 



Physical properties 


. a, 

s dynes/cm 


Values of k for 
atom! zlng-gases 


JP-5 

n-Octane 




TABLE 2. - EFFECT OF PHYSICAL PROPERTIES ON 
RECIPROCAL SMO, AT MACH 1 GASFLOH 


Source 

Relationship 

Gas-phase properties 

Present 

study 

d:' - v'- 33 . p'- 33 . p-°- 7 5, „-0.46 V 2.2S 
32 c K g •'g g ’ m 

Theory 

--1 _ ,.1 .33 0.67 0 u 0 .,0 

D-, V , p . u. , M , V 

32 g K g g m 


Liquid-Jet properties 




Present DZl - p, , y, , o ' 0 ' 44 
study 32 1 ' 

Thpnr v n ' 1 - °- 33 0-67 -0.33 

ineory D „ pj , , a 


Other terms 


Present D'i - D °‘ 2 , g -0 ’ 75 

study 32 0 

Theory D"^ - D°, g° 


Experimental . 
b Calculated. 
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FIGURE 1. - APPARATUS AND AUXILIARY EQUIPMENT. 
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GAS 


INSIDE 
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0.0G4 CM 
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FIGURE 3. - DIAGRAM OF PNEUMATIC TWO-FLUID 
ATOMIZER. 
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LIGHT SCATTER ANGLE, RAD 

FIGURE 4. - TYPICAL SCATTERED-LIGHT 
ENERGY CURVE WITH FOUR CHARACTER- 
ISTIC DIAMETERS. 
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COEFFICIENT, EQ. (3) D n /D- 
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FIGURE 6. - CORRELATION OF SAUTER MEAN DIAMTER, D 32 . 
WITH DIMENSIONLESS GROUP WeRe AND FLUID DENSITY 
RATIO. 
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FIGURE 7. - CORRELATION OF PROPORTIONALITY COEFFICIENT 
k' WITH DIMENSIONLESS GROUP p jV m /g p g . 


10 



(WeRe Pg ^) 0 -' 1 ' 1 <p 1 v r ^/g u g > 0 - 75 

FIGURE 8. - CORRELATION OF NORMALIZED RECIPROCAL SAUTER 
MEAN DIAMETER WITH DIMENSIONLESS GROUPS. 



.1 .2 .3 .<) .6 .8 1.0 2.0 

EXPONENT N n 


FIGURE 9. - CORRELATION OF ROSIN-RAMMLER AND NUKIYAMA 
TANASAWA EXPONENTS N f AND N n , RESPECTIVELY. 


11 


RECIPROCAL SAUTER MEAN DIAMETER 


HELIUM, ATOMIZING-GAS 



FIGURE 10. - COMPARISON OF RECIPROCAL SMD'S FOR WATER 
WITH FUEL-SPRAY CALCULATED VALUES 
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